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Abstract. The coalescence of massive black holes is one of the primary sources 
of gravitational waves (GWs) for LISA. Measurements of the GWs can localize the 
source on the sky to an ellipse with a major axis of a few tens of arcminutes to a few 
degrees, depending on source redshift, and a minor axis which is 2-4 times smaller. 
The distance (and thus an approximate redshift) can be determined to better than 
a percent for the closest sources we consider, although weak lensing degrades 
this performance. It will be of great interest to search this three-dimensional 
"pixel" for an electromagnetic counterpart to the GW event. The presence of 
a counterpart allows unique studies which combine electromagnetic and GW 
information, especially if the counterpart is found prior to final merger of the 
holes. To understand the feasibility of early counterpart detection, we calculate 
the evolution of the GW pixel with time. We find that the greatest improvement 
in pixel size occurs in the final day before merger, when spin precession effects are 
maximal. The source can be localized to within 10 square degrees as early as a 
month before merger at ^ = 1; for higher redshifts, this accuracy is only possible 
in the last few days. 



PACS numbers: 04.80.Nn, 04.30.Db, 04.30.Tv 
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1. Introduction 

One of the most promising sources of gravitational waves (GWs) for LISA is the 
coalescence of massive black hole (MBH) binaries. Observations show that massive 
black holes are present in the center of nearly all galaxies at the present time pQ [5] . 
Hierarchical structure formation teaches that these galaxies formed from the mergers 
of smaller galaxies and their associated dark matter halos [3l[4]. During such mergers, 
dynamical friction and three-body interactions bring the black holes of each galaxy 
together, forming a binary. Eventually the binary is close enough that it decays due 
to gravitational radiation [S]. For 10^ < M < 10^ Mq, these waves lie in the low 
frequency (3 x 10~^ Hz < / < 1 Hz) band of LISA. Merger tree calculations show that 
~ tens of events should be detected by LISA each year [6l |7] . 

The GWs detected by LISA encode a great deal of information about the source 
binary. In general, GWs from MBH binaries are characterized by 17 parameters. 
Since rapid circularization is likely for comparable mass binaries [8] , we consider only 
circular orbits, reducing the parameter set to 15. These parameters are the masses 
of the holes, their spins, the orientation of the binary orbit, the "merger time" and 
"merger phase" (defined by the post-Newtonian approximation; see [9] for details), 
the position of the binary on the sky, and its luminosity distance. LISA estimates of 
the sky position and distance, along with the errors in those measurements, define a 
three-dimensional "pixel" in which the source is expected to be located on the sky. 
It will be interesting to search this pixel for electromagnetic (visible, radio. X-ray) 
counterparts to the GW source. 

It is plausible that no significant electromagnetic activity occurs in conjunction 
with MBH coalescence. In this case, one could imagine searching the pixel for a galaxy 
with a structure that indicates a recent merger. Another possibility is to search for 
a galaxy which has a bulge radial velocity consistent with the GW-measured final 
black hole mass. (This of course assumes that the well-known relation between these 
quantities in the local universe [lOl [11] holds at high redshift and so soon after a 
merger.) 

It is likely, however, that there is some unique EM activity associated with the 
MBH coalescence. The nature of such activity has become a hot research topic in 
the past few years, leading to many possible scenarios [THUS]. For example, if the 
surrounding gas is completely swept away by the binary, there may be no signal 
during the coalescence. Instead, there would be a delayed afterglow when the gas 
later accretes onto the remnant hole [14]. This afterglow may occur years after the 
merger. It is likely, though, that gas won't be totally swept away, leaving enough 
to accrete onto the holes and create variable EM activity during the inspiral phase. 
For example, Armitage & Natarajan showed that for a large mass ratio binary, any 
gas which does remain will be driven in during inspiral, producing an EM signal |15j . 
More recent work by MacFadyen & Milosavljevic shows that periodic variations in the 
Newtonian potential can create a quasi-periodic EM fiux [16]. 

Other scenarios predict transients during or immediately after the coalescence. 
One example relies on the fact that the mass of the black holes is partially radiated 
away by the gravitational waves. This causes a near-instantaneous change in the 
gravitational potential to which the gas will quickly react, producing shocks |17j . 
Another possibility is that the "kick," or momentum imparted to the black hole due 
to an asymmetric emission of GWs [TSl dH HO] , will send the remnant through the 
surrounding gas, again producing shocks [211 [22l [23] . A transient signal might also 



Advanced localization of massive black hole coalescences with LISA 



3 



appear when the GWs are viscously dissipated in the surrounding gas fM\ . 

Finding a counterpart could greatly enhance the science return of MBH 
measurement. For example, counterparts can improve LISA's ability to determine 
certain parameter values. The sky position is correlated with various other parameters, 
particularly luminosity distance and orbit orientation. When it is determined exactly 
by identification of a counterpart, the other parameters can be estimated to greater 
accuracy [25l [26] . Another difficulty with parameter estimation is that the estimated 
masses are actually "redshifted" ; that is, they are equal to the rest frame mass 
multiplied by (1 + The GWs themselves only give the luminosity distance, not the 
redshift, so any decoupling of mass and redshift will be dependent on a cosmological 
model [25]. The counterpart gives the redshift directly. 

Finding the redshift from a counterpart may also allow us to use MBH events 
as cosmological distance measures. Combining the EM-measured redshift with the 
GW-measured luminosity distance creates a Hubble diagram which is calibrated only 
by general relativity [57] [5^ . In reality, such "standard sirens'13 are affected by weak 
lensing of the gravitational waves, which impacts measurements of the luminosity 
distance [28l [29] . Still, the systematics affecting MBH GWs should be completely 
different from those affecting Type la supernova standard candles jSO] [ST] [32| and 
could serve as a useful complement to those sources. 

Counterparts are also useful for studying the astrophysics of the MBH coalescence; 
indeed, the sheer variety of counterpart scenarios shows how uncertain these processes 
are. Specifically, counterparts could give insight into gas dynamics and accretion. For 
instance, GW measurements of the mass and EM measurements of the luminosity 
could be combined to find the Eddington ratio, L/L-^dd [33]. Finally, the counterparts 
could be used to test fundamental physics. If a counterpart features EM variation in 
phase with the gravitational wave signal, the two signals could be compared to test 
the equivalence of photon and graviton propagation speed. Any difference could be 
explained by a nonzero graviton rest mass [13j . 

Given the potential utility of electromagnetic counterparts, it is important to 
understand how well LISA will be able to measure sky position and distance. This 
will give astronomers some indication of LISA capabilities, as well as provide guidance 
regarding LISA mission design decisions. In particular, it is interesting to see how well 
sky position and distance can be measured in advance of merger [34j . The intention is 
to analyze LISA data in real time so that astronomers can point telescopes at the pixel 
to look for and then observe the source. Waiting until after the merger means that 
astronomers will miss any precursor or prompt EM events, which might also make it 
harder to even identify the correct source. 

2. Preliminaries 

To estimate parameter measurement accuracies, we use the Fisher matrix formalism 
[551 [55] ; details for our analysis are given in ^ . The Fisher matrix is defined as 




(1) 



I This term, coined by Storl Phinney and Sean Carroll, is appropriate because GWs can be thought 
of as analogous to sound waves. 
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where h is the GW signal measured by LISA and 6'" are the 15 parameters which 
characterize the signaL The inner product is defined by 



{a\b) = 4 Re 



^ Snif) 



(2) 



Tildes denote Fourier transforms, which we take using the stationary phase 
approximation f361 137] . 5'„(/) is the one-sided noise power spectral density, which 
includes instrument noise [38j and confusion noise from Galactic [39] and extragalactic 
[40j white dwarf binaries. The inverse of the Fisher matrix is the covariance 
matrix, S"^, the diagonal entries of which correspond to (squared) parameter errors. 
The Fisher matrix estimates errors in the "Gaussian approximation," in which the 
parameter posterior probability density is represented by a Gaussian around the most 
likely values. This approximation is known to be good for high signal-to- noise ratio 
|35[ [36], although it is not immediately clear what that means for a signal with 15 
strongly coupled parameters. 

We use only the inspiral portion of the coalescence waveform. This is largely 
because the inspiral can be modeled by the post-Newtonian approximation to general 
relativity [HI \^ offering a convenient parameterization for our analysis. In 
addition, the primary source of information about sky position is the motion of LISA 
around the sun. Therefore, one might guess that the inspiral provides the bulk 
of the localization. However, with the development of numerical relativity, it has 
become possible to perform similar calculations which include the merger phase of 
the waveform. Preliminary results seem to show that the merger could dramatically 
improve sky position and distance determination [44] . 

Our waveforms are carried out to second post-Newtonian (2PN) order in the 
phase, but include only the Newtonian quadrupolc amplitude term. They also include 
the effects of spin-induced precession. Spins precess due to post-Newtonian spin-orbit 
and spin-spin coupling effects. For example, to the order we are considering [451 146j . 
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where nii and TO2 are the black hole masses, M — nii + 7712 is the total mass, 
/u = TO1TO2/M is the reduced mass. Si and S2 are the spin vectors, L is the normal to 
the orbital plane, and r is the orbital separation in harmonic coordinates. On short 
timescales, total angular momentum is conserved, so the orbital plane precesses to 
compensate for the spins. The effect of these precessions on the waveform is to create 
various amplitude and phase modulations. For example, consider the "polarization 
amplitude." It describes the measured amplitude of a GW and is constructed by 
adding the two polarizations, weighted by their antenna pattern functions F'^ and 
F'^, in quadrature [17] : 



A 



poU = — [(1 + (L 



{ty + AiL-nyFtity] 



211/2 



(4) 



where n is the direction to the binary and i G {I, II} denotes one of two detectors 
synthesized from the LISA laser links. Because L precesses, an amplitude modulation 
is imposed on the measured waveform. Other modulations are hidden by the antenna 
pattern functions. 
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Figure 1. "Polarization amplitude" Apoi(t) of the signal measured in detector 
I for a binary with mi = 10® Mq, m2 = 3 X 10^ Mq, and 2 = 1. The curves 
are as follows: solid line, Xl = X2 = 0; dashed line Xi = X2 = 0.1; dotted 
line, XI = X2 = 0.5; and dash-dotted line, Xl = X2 = 0.9. (x = S/m^ is the 
dimensionless spin parameter.) The figure covers the final two years of inspiral. 
Without spin, the only modulation is due to LISA'S orbit around the sun, which 
is the primary source of localization information. Precession adds additional 
modulations, which increase toward merger. 



Figure [T] shows Apoi,i for a two-year cycle. Without precession, the polarization 
amplitude has a period of one year, the time it takes LISA to orbit the sun. This 
orbital modulation (which is hidden in the pattern functions and F^) plays 
a large role in resolving the position of a source on the sky. When precession is 
added, the amplitude is modulated much more strongly. The additional modulation 
is strongest just before merger, when the holes are closer and the post-Newtonian 
precession effects the largest. Similar modulations exist in the relative phase between 
weighted polarizations (the "polarization phase") and the intrinsic post-Newtonian 
phase. Finally, there is an additional precessional correction to the orbital phase |45j . 
All of these effects combine to provide more information about the source, reducing 
parameter errors. Specifically, they break degeneracies between strongly correlated 
parameters, notably the sky position, distance, and orbital orientation. We expect 
the improvements to be greatest when the precession effects are largest, at the end of 
inspiral. 

3. Time evolution of LISA pixel 

The operation of our code is described in [5]. We first select rest frame masses and 
redshift for a binary system. We then randomly choose the spin magnitudes and 
direction, orientation, sky position, and post-Newtonian "merger time" parameter tc- 
For each choice, the code computes the Fisher matrix Tab and the covariance matrix 
Extracting the sky position error ellipse from this requires some care. Defining 
the ellipse such that the probability of the source lying outside it is 63% 
confidence interval) [47j, the major (2a) and minor (26) axes of that ellipse are given 

by 
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the upper sign is for 2a, while the lower sign is for 2b. Here /I at = cos0n, where On 
is the polar angle of the source. 4>n is the azimuthal angle. The bars over the angles 
reflect that they are measured in the barycenter frame, rather than in a frame moving 
with LISA. (See [9] for more details.) 

We find that at merger, the sky position error ellipse has a (median) major axis 
of ~ 15 — 45 arcminutes at z = 1, depending on system mass. This figure increases 
to ~ 3 — 5 degrees at z = 5. The minor axis ranges from ~ 5 — 20 arcminutes at 
z = lto^l — 3 degrees at z = 5. Finally, the luminosity distance errors ADi^/Dj^ 
are - 0.002 - 0.007 at z = 1 and 0.025 - 0.05 at z = 5. It is important to note that 
these luminosity distance errors do not include the effects of weak lensing. Lensing 
due to intervening matter will magnify or demagnify the waves, causing an incorrect 
measurement of luminosity distance to the few percent level [28l [29| . 

We take an area of 10 square degrees as a useful benchmark for a "good" 
localization. This is roughly the field of view of the Large Synoptic Survey Telescope 
(LSST) [48j. For the luminosity distance, we consider a source to be well-localized 
when the distance error calculated from only the GWs drops below the error due to 
weak lensing, which we can roughly approximate as 5% for this purpose. Then for 
low redshift (z = 1), sources are well-localized as far back as a month before merger 
for most masses. At z = 3, the sky position ellipse is smaller than 10 square degrees 
a few days before merger in the smallest two-thirds of the masses we consider. The 
luminosity distance error drops below the weak lensing limit a few days to a week 
before merger. Finally, at z = 5, the error ellipse is the size of the LSST field or 
smaller at most a day before merger, and only then in a few cases. The luminosity 
distance errors are worse than the lensing limit at all times before merger. 

Figure [2] shows the evolution of the major axis and minor axis of the sky position 
error ellipse for mi = IO^Mq, m2 = 10^ Mq, and z = 1. Each histogram represents 10** 
different choices of the random parameters. As the time until final merger decreases 
from 1 month to 1 day, the median error gradually decreases. In addition, the shape 
of the distribution changes slightly, becoming less peaked. However, a drastic change 
occurs between 1 day and days (end of inspiral). In just this day, the median is 
reduced by half an order of magnitude, and the shape of the distribution changes 
drastically. For the major axis, the distribution becomes almost fiat over a certain 
range and displays a slightly bimodal structure. 

Figure [3] shows the median values of our distributions for different masses. We see 
a trend holding over all masses: gradual improvement with time until a sharp change 
the last day. The change is greatest for the lowest mass systems. 

These results agree quite well with those of Kocsis et al [3H] , except their study 
does not show a sharp improvement in localization during the final day of inspiral. It 
turns out that this eS'ect is due mostly to spin precession, which Kocsis et al do not 
include. This late impact of precession is to be expected following our discussion of 
Figure [Hand is confirmed in Figure ID By breaking correlations between parameters, 
precession gives a factor of 2 — 7 improvement in final parameter estimation errors. 
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Figure 2. Distribution of the major axis 2a (top) and minor axis 2b (bottom) of 
the sky position error ellipse for lO** binaries with mi = 10® Mq, m2 = 10^ Mq, 
and 2 = 1 at different times before merger. Reading from left to right, the number 
of days until merger is (solid line), 1 (dashed line), 2 (dash-dotted line), 4 (dotted 
line), 7 (solid line again), 14 (dashed line), 21 (dash-dotted line), and 28 (dotted 
line). The largest change occurs in the final day before merger. 



However, since most of this improvement takes place in the final days before merger, 
precession does not help very much with advanced localization. This effect may, 
however, impact the frequency of data downlink between LISA and the ground: As 
a source approaches merger, it may be useful to continuously downlink data to take 
advantage of the rapid improvement in localization. 

As a final study, we also investigate the position dependence of the localization 
pixel. Rather than randomly choosing the full sky position, we randomly choose just 
one of the polar (^jv) and azimuthal {(f>N) angles, while varying the other in a controlled 
fashion. It turns out that only the polar dependence is meaningful due to a particular 
feature of our formalism: We randomly choose the final merger time within LISA's 
mission window, meaning that even with a specifically chosen ^j^, the relative azimuth 
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Figure 3. Median values of 2a, the sky position ellipse major axis, as a function 
of time before merger. Each point is taken from a Monte Carlo sample of 10* 
binaries with fixed masses (given in units of Mq) and z = 1. Data were only 
output at the marked points; the lines are there to guide the eye. 



between LISA and the source at merger, S(j) — (pM — 4''LisK{ic), is random. This mostly 
washes out the azimuthal dependence. Other authors HOI US] niake a different 
choice, in which the position of LISA at merger is fixed. This allows them to study 
the position dependence of localization at a particular time of the year; consequently, 
their results for both polar and azimuthal angle dependence are different from ours. 

We find that the sky position axes and area are maximal for sources in the ecliptic 
plane and minimal near the poles, at cos^^v — ±0.8. The distance has a three-peaked 
structure, with one peak at the ecliptic plane and the others very close to where the 
sky position errors are minimized. Converting to Galactic coordinates, we find that 
some of the best localized sources lie outside of the Galactic plane. This suggests that 
the most easily localized binaries will not tend to lie in regions of the sky with thick 
foregrounds or high levels of extinction. 
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Figure 4. Medians of 2a (solid lines) and 2b (dashed lines) as a function of time 
for binaries with mi = 3 X IO^Mq, m2 = 10^ Mq, and 2 = 1. Spin-induced 
precession was included in the waveform for the crosses and neglected for the 
circles. 



4. Summary and future work 

Observing electromagnetic counterparts to massive black hole binary coalescences can 
be very useful in the contexts of parameter estimation, cosmology, astrophysics, and 
fundamental physics. To observe a counterpart, and to get the most science out of it, 
it is useful to be able to localize the source prior to merger. Our study shows that 
advanced localization of sources with LISA should be possible at low redshift {z = 1), 
but begins to be difficult at higher z. The effects of precession break degeneracies 
between parameters, dramatically improving estimates in the final days before merger. 
Finally, at least some of the best located sources will be out of the Galactic plane, 
improving our chances of identifying counterparts. 

These results are not the final story on LISA's ability to measure and 
localize MBH sources. For example, it is useful to test whether the Gaussian 
approximation/Fisher matrix approach is actually valid for a signal with 15 
parameters. This problem is currently being approached both by analytic analysis 
|51| and by comparing the Fisher results to errors generated by a Markov Chain 
Monte Carlo exploration of the parameter posterior distribution function [35] ■ Another 
approximation worth checking is the stationary phase method for obtaining Fourier 
transforms; this is known to be good for signals without precession [5 3) but seems 
to smooth out sharp features caused by strong precession [54]. Finally, there are 
various extra pieces which can still be added to the waveform, including higher 
post-Newtonian phase terms, higher order precession equations [55l [56] , higher post- 
Newtonian amplitude terms, and harmonics beyond the quadrupole. The latter two 
terms, often called simply "higher harmonics," have already been shown to decrease 
parameter errors in the absence of spin precession [57|, I50|, 158] . The improvement 
seems to be somewhat orthogonal to the improvement caused by precession, but a 
more complete investigation, including time dependence, is necessary. We and our 
collaborators are currently investigating these issues. 



Advanced localization of massive black hole coalescences with LISA 



10 



Acknowledgments 

Our work on LISA measurement analysis has been supported by NASA Grants 
NNG05G105G and NNX08AL42G, as well as NASA Contract No. 1291617 and the 
MIT Class of 1956 Career Development Fund. We also gratefully thank the organizers 
of the Symposium for the travel support provided to RNL and for the outstanding 
meeting that they put together. 

References 

[1] Kormendy J and Richstone D 1995 Ann. Rev. Astron. Astrophys. 33 581 

[2] Magorrian J et al 1998 Astron. J. 115 2285 

[3] Madau P 2007 Preprint arXiv:0706.0123 

[4] Di Matteo T, Colberg J, Springel V, Hernquist L and Sijacki D 2008 Astrophys. J. 676 33 

[5] Begelman M C, Blandford R D and Rees M J 1980 Nature 287 307 

[6] Sesana A, Volonteri M and Haardt F 2007 Mon. Not. R. Astron. Soc. 377 1711 

[7] Micic M, HoUey-Bockelmann K, Sigurdsson S and Abel T 2007 Mon. Not. R. Astron. Soc. 380 
1533 

[8] Peters P C 1964 Phys. Rev. 136 B1224 

[9] Lang R N and Hughes S A 2006 Phys. Rev. D 74 122001 

[10] Ferrarese L and Merritt D 2000 Astrophys. J. Lett. 539 L9 

[11] Gebhardt K et al 2000 Astrophys. J. Lett. 539 L13 

[12] Dotti M, Salvaterra R, Sesana A, Colpi M and Haardt F 2006 Mon. Not. R. Astron. Soc. 372 
869 

[13] Kocsis B, Haiman Z and Menou K 2008 Astrophys. .J. 684 870 

[14] Milosavljevic M and Phinney E S 2005 Astrophys. J. Lett. 622 L93 

[15] Armitage P J and Natarajan P 2002 Astrophys. J. Lett. 567 L9 

[16] MacFadyen A I and Milosavljevic M 2008 Astrophys. J. 672 83 

[17] Bode N and Phinney E S 2007 APS Abstract http://meetings.aps.org/link/BAPS.2007.APR. Sl.10l 

[18] Fitchett M J 1983 Mon. Not. R. Astron. Soc. 203 1049 

[19] CampaneUi M, Lousto C O, Zlochower Y and Merritt D 2007 Phys. Rev. Lett. 98 231102 

[20] Schnittman J D, Buonanno A, van Meter J R, Baker J G, Boggs W D, Centrella J, Kelly B J 

and McWilliams S T 2008 Phys. Rev. D 77 044031 

[21] Lippai Z, Frei Z and Haiman Z 2008 Astrophys. J. Lett. 676 L5 

[22] Shields G A and Bonning E W 2008 Astrophys. J. 682 758 

[23] Schnittman J D and Krolik J H 2008 Astrophys. J. 684 835 

[24] Kocsis B and Loeb A 2008 Phys. Rev. Lett. 101 041101 

[25] Hughes S A 2002 Mon. Not. R. Astron. Soc. 331 805 

[26] Holz D E and Hughes S A 2005 Astrophys. J. 629 15 

[27] Schutz B F 1986 Nature 323 310 

[28] Wang Y, Holz D E and Munshi D 2002 Astrophys. J. Lett. 572 L15 

[29] Dalai N, Holz D E, Chen X and Frieman J A 2003 Astrophys. J. Lett. 585 Lll 

[30] Phillips M M 1993 Astrophys. J. Lett. 413 L105 

[31] Riess A G, Press W H and Kirshner R P 1995 Astrophys. J. Lett. 438 L17 

[32] Wang L, Goldhaber G, Aldering G and Perlmutter S 2003 Astrophys. J. 590 994 

[33] Kocsis B, Frei Z, Haiman Z and Menou K 2006 Astrophys. J. 637 27 

[34] Lang R N and Hughes S A 2008 Astrophys. J. 677 1184 

[35] Finn L S 1992 Phys. Rev. D 46 5236 

[36] Cutler C and Flanagan E E 1994 Phys. Rev. D 49 2658 

[37] Poisson E and Will C M 1995 Phys. Rev. D 52 848 

[38] Larson S L, Hiscock W A and Hellings R W 2000 Phys. Rev. D 62 062001 

[39] Nelemans G, Yungelson L R and Portegies Zwart S F 2001 Astron. Astrophys. 375 890 

[40] Farmer A J and Phinney E S 2003 Mon. Not. R. Astron. Soc. 346 1197 

[41] Blanchet L 2006 Liv. Rev. Rel. 9 4 

[42] Blanchet L, Damour T, Iyer B R, Will C M and Wiseman A G 1995 Phys. Rev. Lett. 74 3515 

[43] Will C M and Wiseman A G 1996 Phys. Rev. D 54 4813 

[44] Babak S, Hannam M, Husa S and Schutz B 2008 Preprint ar Xiv:0806.1591l 

[45] Apostolatos T A, Cutler C, Sussman G J and Thorne K S 1994 Phys. Rev. D 49 6274 

[46] Kidder L E 1995 Phys. Rev. D 52 821 



Advanced localization of massive black hole coalescences with LISA 



11 



[47] Cutler C 1998 Phys. Rev. D 57 7089 

[48] Tyson J A et al 2002 Proc. SPIE 4836 10 

[49] Kocsis B, Haiman Z, Menou K and Prei Z 2007 Phys. Rev. D 76 022003 

[50] Trias M and Sintes A M 2008 Phys. Rev. D 77 024030 

[51] Vallisneri M 2008 Phys. Rev. D 77 042001 

[52] Cornish N J and Crowdcr J 2005 Phys. Rev. D 72 043005 

[53] Droz S, Knapp D J, Poisson E and Owen B J 1999 Phys. Rev. D 59 124016 

[54] Cornish N J, O'Sullivan S, Lang R N and Hughes S A, in preparation 

[55] Fayc G, Blanchct L and Buonanno A 2006 Phys. Rev. D 74 104033 

[56] Blanchct L, Buonanno A and Fayc G 2006 Phys. Rev. D 74 104034 

[57] Arun K G, Iyer B R, Sathyaprakash B S, Sinha S and Van Den Broeck C 2007 Phys. Rev. D 76 
104016 

[58] Porter E K and Cornish N J 2008 Phys. Rev. D 78 064005 



